Rational Orbits around Charged Black Holes 
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We show that all eccentric timelike orbits in Reissner-Nordstrom spacetime can be classified 
using a taxonomy that draws upon an isomorphism between periodic orbits and the set of rational 
numbers. By virtue of the fact that the rationals are dense, the taxonomy can be used to approximate 
aperiodic orbits with periodic orbits. This may help reduce computational overhead for calculations 
in gravitational wave astronomy. Our dynamical systems approach enables us to study orbits for 
both charged and uncharged particles in spite of the fact that charged particle orbits around a 
charged black hole do not admit a simple one-dimensional effective potential description. Finally, 
we show that comparing periodic orbits in the RN and Schwarzschild geometries enables us to 
distinguish charged and uncharged spacetimes by looking only at the orbital dynamics. 



I. INTRODUCTION 

Black holes are useful large-scale laboratories for test- 
ing general relativity in strong gravitational fields. While 
directly observing black holes proves difficult because 
they emit no electromagnetic radiation, black hole pairs 
can be detected via the gravitational radiation they may 
emit. The terrestrial network of interferometric gravita- 
tional wave detectors and the proposed space-based de- 
tector LISA (Laser Interferometer Space Antenna), are 
expected to detect gravitational radiation and launch an 
era of gravitational wave astronomy that would make use 
of direct observations of black holes. 

Because gravitational waves are shaped by the mo- 
tion of massive celestial objects, extracting astrophysi- 
cally meaningful information from them requires a com- 
prehensive theoretical understanding of the sources' un- 
derlying dynamics 16 . For example, one popular pro- 
cessing method, matched filtering, makes use of a tem- 
plate signal, generated using theoretical predictions, to 
find signals in noisy detector output [7J [5]. Template 
signal generation is an example of the type of compu- 
tationally expensive process encountered in gravitational 
wave astronomy when studying aperiodic orbits in the 
strong-field regime. 

A method for approximating aperiodic orbits with pe- 
riodic orbits, which could cut down significantly on com- 
putational expenses, was introduced in an earlier paper 
[9] . The approximation method takes the form of a tax- 
onomy that assigns to each periodic orbit a rational num- 
ber. By virtue of the fact that the rationals are dense, the 
taxonomy can be used to approximate aperiodic orbits 
with periodic orbits to arbitrary precision. Because peri- 
odic orbits might have Fourier series that converge more 
rapidly than those of aperiodic orbits, and because for 
periodic orbits the evolution of a geometry's conserved 
quantities may be interdependent, calculations pertain- 
ing to periodic orbits might be less computationally in- 
tensive than those for aperiodic orbits [5]. 

The taxonomy was applied to the Kerr geometry in [9 



and [TU] and to black hole pairs in [IT] and [12] • Here, 
we will extend the approach to the Reissner-Nordstrom 
(RN) solution to the Einstein field equations, which de- 
scribes the gravitational field of a static, non-rotating, 
electrically charged, spherically symmetric body [T3j . 
Studying a geometry of this type is less astrophysically 
motivated than studying its electrically neutral counter- 
part. Were it to form in spite of the fact that the electro- 
magnetic repulsion in compressing an electrically charged 
mass is greater by about 40 orders of magnitude than 
the gravitational attraction, it would neutralize its own 
electric charge if enough opposite charge were available. 
Nonetheless, in the spirit of being prepared for the un- 
expected, and in support of the ambitious gravitational 
wave experiments coming online, we will not presump- 
tively exclude any possible sources. Were gravitational 
waves from a charged black hole candidate detected ex- 
perimentally, a thorough understanding of the RN orbits 
could help identify their source. 



At the other extreme, microscopic black holes that 
might form in accelerator experiments may be charged. 
A pair of black holes that scatter and then evaporate 
might be described by a scattering amplitude that is a 
sum over these classical paths. The solutions we describe 
might find application in particle physics as well as as- 
trophysics. 



Both Kerr and RN orbits bear many qualitative sim- 
ilarities to orbits in Schwarzschild spacetime, including 
"zoom- whirl" behavior [7J [2], in which the test parti- 
cle zooms away from the central mass quasielliptically to 
successive apastra separated by nearly circular whirls. It 
is this behavior which the taxonomy exploits. As was the 
case for Kerr orbits before this taxonomy was introduced, 
no unifying framework for making general claims about 
orbits has been applied to the RN spacetime before. 



II. REISSNER-NORDSTROM EFFECTIVE 
POTENTIAL 

We use the effective potential formulation of the 
RN spacetime to calculate various interesting dynami- 
cal properties of the geometry. In an appendix we detail 
the Hamiltonian formulation that is used to generate the 
orbits pictured throughout the paper. We begin with the 
RN metric, 
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where the horizon function A is 
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We have assumed that the central magnetic charge is 
zero, which would otherwise change the horizon func- 
tion's third term to (Q 2 + P 2 )/r 2 . We use geometrized 
units and measure r and the central charge Q in units of 
M. 

A first constant of motion is always 



flf,, 



dx* dx v _ 
dr dr 

for timelike geodesies 
for null geodesies 



Explicitly, for timelike orbits, this gives 
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r 2 + r 2 9 2 + r 2 sin 2 9 ip 2 = -1, (4) 



where an overdot indicates differentiation with respect to 
an affine parameter r. 

To derive expressions for the conserved quantities we 
write the Lagrangian for a charged particle, in which the 
final term is derived using the fact that the only nonva- 
nishing component of the vector potential A^ is Aq [TS] : 



£ = i(-AP. 
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-r 2 9 2 + r 2 sin 2 i 
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Q* is the charge per unit mass of the test particle. 
Geodesies in the RN geometry are constrained to a plane. 
Since the geometry is spherically symmetric, we can set 
9 = tt/2 and 9 = so that every geodesic is equatorial. 

Using 



d_ fdC\ _d£ 
dr \dq J dq 

yields the constants of motion: 



p t — -At = -E 
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P v = r 2 ip = L. (7) 



Applying Q to Q and simplifying gives us 

1 „ 

-f + V cS = £ eS , 

where 
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Due to the geometry's spherical symmetry, the ef- 
fective potential is one-dimensional - just as in the 
Schwarzschild case - but with one caveat: the test par- 
ticle's charge Q* is coupled to the energy E. This pro- 
hibits us from writing V e s in a form that is independent 
of energy without also writing £ c ff in a form that is not 
constant. First we will consider the case when Q* = so 
that we can define a true one-dimensional effective poten- 
tial. Applying this condition to (foj) gives us the effective 



RN potential for an uncharged particle: 
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We return to the charged test particle in { IV 



III. UNCHARGED PARTICLE ORBITS IN RN 
SPACETIME 

A. Bounds on Q 

For an uncharged particle, the RN geometry is quali- 
tatively similar outside the horizon to the Schwarzschild 
geometry [T5] . But because the RN V c g is larger than the 

Schwarzschild V c g by a factor of ^ + w 2ri , a thorough 
quantitative analysis is necessary to understand the or- 
bital dynamics. Our goal in this section is to determine 
the bounds on L and Q that yield periodic orbits. 

Determining bounds on Q is simple: owing to the RN 
geometry's peculiar horizon structure, only certain values 
of Q are realistic [T3]. The null hypersurfaces are given 
by 



A(r) = 1 - 



= 0, 



which when solved for r gives us the horizons: 



r± = 1 ± y/1 - Q 2 . 



(12) 



(13) 



When Q 2 > 1, the geometry is a naked singularity. 
The condition Q 2 = 1 describes an extremal RN black 
hole, which is highly unstable due to the fact that adding 
any mass at all makes it an undercritically charged black 
hole. This leaves the undercritically charged geometry 
(Q 2 < 1) as the most realistic scenario. Because Q 
only appears in the potential as Q 2 , the undercritically 
charged case is equivalently given by — 1 < Q < 1. 

Q affects the shape of the potential as shown in Fig- 
ure [I] Note that the Q — case (solid) is simply the 
Schwarzschild effective potential. When comparing the 
RN potential to the Schwarzschild potential, we see that 
the factor of +Q 2 /2 in the r~ 2 term is the reason for 
the heightened peak and that the r term makes the 
potential blow up at zero. 
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The Schwarzschild V c g has at most two extrema. Be- 
cause the RN Vqs is a quartic, it can have three extrema. 
This opens up the possibility of three circular orbits - 
two stable and one unstable - outside the horizon. The 
existence of two stable circular orbits would imply that 
there are two regions in which we can find bounded or- 
bits. There would be two stable circular orbits if both 
minima were to lie outside the outer horizon, r+. Fig- 
ure [2] depicts the Schwarzschild and RN potentials for 
L = 3.2. The dashed potential is the Schwarzschild V e g 
and the Schwarzschild horizon is given by the dashed 
vertical line (r = 2). The RN V e g and external horizon 
are solid. For this choice of L and Q there are only two 
extrema outside the RN horizon, one stable and one un- 
stable. We want to determine if there are any L and Q 
for which the external horizon is closer to the singularity 
than the inner minimum of the potential. This would 
yield two stable circular orbits. 
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FIG. 2: The RN and Schwarzschild potentials with horizons 
for each and a non-physical internal horizon (dotted) that rep- 
resents a scenario in which there are multiple stable circular 
orbits 



Our first step is to find the radii r c at which the effec- 
tive potential has stationary points, given by 



dVeg = 1_ _ Q 2 +L 2 + 31? _ 2Q 2 L 2 = Q 

yt Zi iy**j rv*-± fr*U 



dr 



We rewrite this as 



(Q 2 + L 2 )r 2 c +3L 2 r c -2Q 2 L 2 = (15) 



FIG. 1: Kff for Q = 0, 0.4, 0.8, and 1, from bottom to top and take itg d i scrirnmant to gct 
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When D > 0, 



has 3 distinct, real roots, and Vvfr has 



d" f * nas o uisuiici, real roots, aiiu v c ff 
3 extrema. If the condition D > is ever true when the r c 
are all > r + , the external horizon, we will have two stable 
circular orbits. The solutions of Equation (151 for r c do 



not provide much insight, so we will not reproduce them 
here. Instead we represent the solutions graphically. 

The disjointedness of the regions in parameter space 
that satisfy each of the above conditions is demonstrated 
in Figure [31 The regions with horizontal hatching are 
those in which D > and the regions with vertical hatch- 
ing are where the smallest circular orbit r c > r + . Note 
that the vertically hatched region does not presume the 
existence of three extrema. There is no overlap between 
the two regions when — 1 < Q < 1 and Q* = 0, so in 
these conditions we never see three extrema outside the 
event horizon. Therefore, there is always at most one sta- 
ble circular orbit. This clarifies the region of the effective 
potential in which we find periodic orbits and confirms 
that for Q* = 0, there is always at most one stable cir- 
cular orbit. 
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FIG. 3: Regions in which the potential has three extrema 
(horizontal hatching) and those in which the smallest ex- 
tremum is outside the horizon, (vertical hatching) 
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FIG. 4: V a e for L = 3.1, 3.4, 3.8, and 4.0, from bottom to 
top, with Q fixed. 



where ISCO stands for "Innermost Stable Circular Or- 
bit" and IBCO for "Innermost Bound Circular Orbit". 
Lisco is the- lowest value of L for which the potential has 
a local minimum, and therefore marks the first appear- 
ance of a stable circular orbit. For L < LiscOj all orbits 
will plunge into the black hole, so Xisco sets the lower 
limit on bound orbits. Libco marks the first appearance 
of an unstable circular orbit that is energetically bound. 
It sets the upper limit only in the sense that we expect 
to see the most zoom-whirl behavior in the strong-field 
when such an unstable bound orbit comes into play [S]. 
Orbits will whirl more as they roll up the potential to- 
wards the unstable bound circular orbit. 

To derive expressions for Lisco and Xibco we start 
with the conditions for circular orbits. Equation (fSj) , 
written with V c s as it appears in Equation (11), is 
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(18) 



B. Bounds on L 

We have established that in the undercritically-charged 
RN geometry for a neutral particle, for any L there is at 
most one stable circular orbit. Figure [4] depicts the po- 
tential for various L with Q fixed. We will now determine 
the bounds in the undercritically-charged RN geometry 
that specify the region in which we find most zoom-whirl 
behavior [9j, 



-kisco < L < LibcO: 



(17) 



The two conditions for circular orbits are r = and r = 0, 
which imply that 



V cS = 
dV eS 



E 2 -l 



dr 



= 0. 



The first condition, when solved for L 2 , gives us 



L 2 ={^-l\r 2 



(19) 



(20) 



We may rewrite j¥- as 
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Applying the conditions for circular orbits in Equation 
( 19 1 gives us 



dVes = f-L 2 A' (E 2 

dr { r 3 A 2 I 2 
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Solving for i 2 gives us 
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Equating Equations (20 1 and (25) and solving for E 2 
gives 



E 2 c (r c ,Q) = 



2A 2 



2A-rA r 



(26) 



Putting Equation ( 26 ) into Equation ( 20 1 and solving for 
L 2 yields 



L 2 c (r c ,Q)=r 2 



rA' 



2A - r A' 



(27) 



With L c and I? c in hand, each of which is depicted in 
Figure [5j our next task is to calculate nsco an d nsco- 
Because at infinity the potential is 0, tibco is simply the 
radius of unstable circular orbit r c such that T4ff(r c ) = 0. 
We first solve 
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for L 2 to obtain 
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(29) 



Setting this equal to L 2 gives us an expression we may 
use to find riBCO : 
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Solving 


this for r 


gives us 
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G(Q) = v/l28 - 144Q 2 + 27Q 4 + 3V-96Q 6 + 81Q 



imco is then L c with hbco m place of r: 
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We determine nsco by taking V£g = (differentiated 
with respect to r), or 



d 2 V cS 3L 2 A 2L 2 A' (L 2 
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Solving for L 2 gives us 



r,2 



J V"=0 4AV _ g A _ A „ r 2 

Setting iy//_ equal to L 2 and solving for r gives us 
4-3Q 2 



(37) 



nsco = 2 



H(Q) 



H(Q) 



(38) 



H{Q) = y/8 - 9Q 2 + 2Q 4 + V5Q 4 - 9Q 6 + 4Q», 
which when plugged into L c gives us 

~Q r ISCO + r ISCO 



r 2 
-°isco 



2Q 2 - 3r ISCO + rf sco ' 



(39) 



Hbco and nsco are plotted in Figure [6j 

We can confirm that the presence of three station- 
ary points in the potential is tied to Liqqo by check- 
ing whether the regions in parameter space D > and 
L 2 > L 2 sco coincide. Figure [H shows that this is the case 
within our bounds for Q. 

Having established bounds on Q and L for the Q* = 
case, we can taxonomize all zoom-whirl behavior in the 
strong field bounded by £isco < L < LibcO: as we show 
in SJV] First, we briefly discuss homoclinic orbits and then 
the Q* ^ case. 



C. Homoclinic Orbits 

For any L such that Lisco < L < £ibcOi there is 
an unstable circular orbit called a homoclinic orbit. The 
homoclinic orbit is the separatrix between orbits that 
plunge to the horizon and those that do not [16] . It pro- 
vides us with the infinite-whirl limit and is therefore an 
important landmark in the orbital landscape. Because 
homoclinic orbits are central to strong-field dynamics 
O [12j [TBI Hi]) we include the orbital plot of a homo- 
clinic orbit for an uncharged particle in RN spacetime in 
Figure [8] 
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FIG. 5: L 2 and E\ as functions of r for Q = 1, Q — 0.5, and 

Q = o. 
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FIG. 7: Regions in which L 2 > L 2 sco (horizontal hatching) 
and where the discriminant D > (vertical hatching) . When 
we impose the bound — 1 < Q < 1 the two regions are equiv- 
alent, which demonstrates that only when L 2 > L 2 SC o does 
the potential have three extrema. 



15; 

10 ; 
5; 

o; 

-5; 
10; 

-15 : 



L=3.3,Q=1.0,Q*=0 




q=oo 



E = 0.994420 

-15 -10 -5 5 10 15 



FIG. 8: A homoclinic orbit for Q = 1.0, Q* = 0, L = 3.3 
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FIG. 6: hbco (solid) and nsco (dotted) as a function of Q 



IV. CHARGED PARTICLE ORBITS IN RN 
SPACETIME 



Having discussed the uncharged test particle scenario, 
we will now move on to orbits in which the particle has 
nonzero charge. While the Q* 7^ effective potential 
(Equation Q) depends on both L and E, it still only 
depends on one dynamical variable, r. Because there is 
no dependence on (9, ip or t, every bound orbit has fixed 
apastra and periastra that are functions of Q, L, and 
E, so we still find periodicity. As such, we can make 
use of this potential to study orbits for charged particles. 
However, the method for finding periodic orbits is more 
complicated than in the Q* = case, because given a 
particular effective potential we may no longer choose £ c s 
so as to produce a periodic orbit, because the potential 
is coupled to E. 
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FIG. 9: Veft for Q, = {0,0.1,0.2,0.3,1.0} from bottom to 
top, with £ sS = ~ 1 + g2 = 0.99, Q = 1, L = 3.5 



A. Identifying Regions with Bounded Orbits 

Because we are considering only test particles we as- 
sume that the particle mass m is << M, which implies 
that Q* < (m/M) < 1. As before, we first want to deter- 
mine whether the potential ever has three extrema out- 
side the horizon, which we do by determining if it is ever 



dr 



IS 



simultaneously the case that the discriminant of 
positive - which means that there are three extrema 
and that the smallest of the extrema is at a radius > r + . 
The region of parameter space in which each condition is 
true is depicted in Figures [TO] and [TTJ respectively. The 
discriminant here is not the same as in Equation (16 1, 
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FIG. 10: Regions of positive discriminant for L 
{3.2,3.5,3.8}; hatching: (diagonal, vertical, horizontal); E 
0.98 



but rather 
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Each region in Figure [TT] is a subset of the correspond- 
ing region in Figure [lOj which means that for these values 
of L, there are in fact regions in which the potential has 
three extrema - two minima - outside the horizon. Note 
that this only occurs for near-extremal geometries and 
only when Q and Q* have like sign, but that for all Q» 
there exists a Q for which there are three external ex- 
trema. 

This is not the complete picture. The existence of two 
minima outside the horizon does not imply that there are 
multiple stable circular orbits, because unlike in the Q* = 
case, £ c ff is fixed by the E we chose for plotting the 
potential. For there to be two stable circular orbits such 
that an energetic particle could "roll" from one stable 
circular orbit over the local maximum and into the stable 
circular orbit, the following condition must be satisfied: 



V e s(ri,L,Q,q,E) = V r ff(r 2 ,L,Q,g, E), 



(41) 



where r\ and ri are the radii of stable circular orbits. 

This is a special case of when the effective potential 
has two minima and £ e g cuts through the more central 
of them in such a way that the periastron of the orbit it 
defines is outside the horizon. If the general case never 
occurs, the special case is impossible as well. We check 
for the general case using a parameter plot, shown in Fig- 
ure 12 Were the three regions shown in the plot ever to 



overlap, we would have parameters for which the poten- 
tial has two minima such that £ e s cuts across the first 
minimum to yield periodic orbits. The regions in this 




FIG. 11: Regions in which the potential has three extrema 
that are all > r+, for L = {3.2,3.5,3.8}; hatching: (diagonal, 
vertical, horizontal); E = 0.98 



plot do not overlap for Q = 1.0, L = 3.6. Regenerating 
the plot for various Q and L reveals the effect of changing 
each on the shape of each region and makes it clear that 
the three regions never overlap. We cannot expect to 
find any choice of parameters that yields a second region 
of bounded orbits in the potential's more central mini- 
mum. This outcome does not preclude the existence of 
bounded orbits in the second minimum, where there arc 
no obstacles to applying the taxonomy. 



B. Bounds on L 

The final step in characterizing the charged particle 
RN geometry is to calculate Lisco an d ^ibco- We first 
write V^' ff and V£g. 



_ 2EQ,Q-2 SL 2 + 3Q 2 -3Q,Q 2 
WL 2 Q 2 
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The condition for circular orbits is V£ s — 0; solving this 
for L 2 and equating it with the result of solving V e s = 
for L 2 yields an expression we may solve for E 2 



given m 




FIG. 12: Q — 1.0; L — 3.6. Each point in the parameter plot 
defines an effective potential and £ e a- The regions shown 
are those where the resulting orbits' periastra are at a radius 
outside the horizon (horizontal hatching), where the perias- 
tra are inside the potential's smaller minimimum (vertical 
hatching), and where the discriminant is positive (diagonal 
hatching). The coincidence of the three regions gives us the 
subset of parameter space in which we can find bounded or- 
bits in the first of the potential's two minima 
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Equation (44). Plugging E 2 
yields L 2 , also given below. 



into Ly, =0 and simplifying 



E 2 e (r,Q,Q* 



L 2 c (r,Q,Q«) = 
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(44) 



The ISCO exists at the inflexion point, given by V"(r) = 
0. So the solution of L v , =0 = L v „ =0 for r gives us nsco- 
Likewise, the solution of L v , =0 = L v=0 for r gives us 
?"ibco Plugging each of these into L\ then yields Libco 
and Lisco- Because these solutions do not provide any 
additional insight, we do not reproduce them here. 

Having defined the bounds in which we find periodic 
orbits, we can move on to applying the taxonomy to or- 
bits in RN spacetime. The taxonomy can be used to 
compare the parameters that yield orbits of a given q in 
the Q 7^ 0, Q* = 0; Q ^ 0, Q* ^ 0; and Q = cases, as 
we will now show. 



V. PERIODIC TABLES OF ORBITS IN RN 
SPACETIME 

A. Overview of the Taxonomy 

Before discussing the taxonomy as it applies to the RN 
spacetime, we will summarize its salient features, which 
were presented in an earlier paper [5] . As mentioned, or- 
bits to which we can apply this taxonomy appear within 
the Libco an d Lisco of the geometry, which define the 
region of parameter space in which the potential accom- 
modates bounded orbits. 

The taxonomy assigns to each orbit a distinct ratio- 
nal number using a scheme that takes advantage of the 
fact that in the strong-field regime, every periodic orbit 
exhibits certain clearly visible topological characteristics. 
In this section we will discuss how to assign to each or- 
bit a rational number based on the orbit's topological 
features. 

Each periodic orbit is associated with a rational num- 
ber g, defined as 



q = w - 



(45) 



where w, v, and z are integers. Each of these integers 
corresponds to a specific topological characteristic of a 
given periodic orbit. The most easily visualized is "z," 
the number of leaves, or "zooms" in the particle's orbit. 
In its path around an RN black hole, our test particle 
will trace out a number of leaves before closing. Figure 
13 depicts orbits with various z values. 

The integer w defines the number of "whirls" the par- 
ticle makes in its path from apastron to periastron to 
the subsequent apastron. To understand this, note that 
every object travels at least a full 2ir around the central 
black hole. The number of whirls is defined as the ad- 
ditional integer number of 2-7T executed beyond this. In 
other words, the number of extra turns around the center 
of the geometry gives us the value of w. Figure [14] shows 
orbits with various w values. 

We require a third number v, the "vertex" number, 
to distinguish between orbits that have equal z and w 
but are geometrically different nonetheless. This is easily 
seen in Figure \TE\ in which both orbits have z = 4, w = 1, 



but where we see that the particle can skip leaves in its 
motion from apastron to apastron. 

We label successive apastra of a periodic orbit with in- 
tegers, counting the starting apastron of the orbit as 
and increasing in the same rotational sense as the orbit 
(counterclockwise for prograde orbits, clockwise for ret- 
rograde orbits), as shown in Figure 15 In general, any 



periodic orbit with z > 2 can skip any number of vertices 
less than z when moving between apastra. We define v 
to be the index of the first vertex hit by the orbit after 
v = (the bounds on v are therefore 1 < v < z — 1). 
When z = 1, we define v = 0, which is the only sensible 
choice for v because it implies that that successive apas- 
tra for single-leaf (z = 1) orbits are actually the same 



single apastron (see Figure 13) 



Finally, we must address the degeneracy that arises 
when the quotient v/z is a reducible fraction. For a given 
w, there are multiple choices of z and v that describe 
the same orbit; for instance, when we have q = lj the 
particle skips every other apastron and never hits vertices 
1 or 3, which means it closes after only two leaves. This 
is equivalent to a q — l| orbit. To avoid this issue, we 
require that v and z be relatively prime. The bounds on 
v are therefore 



1 < v < z — 1 if z > 1 and z, v are relatively prime 
v = if z = 1. (46) 

The rational q defines all of the topological features of any 
closed equatorial orbit and corresponds to the precession 
of the orbit beyond a Keplerian ellipse. 

The utility of the taxonomy becomes apparent when 
we compare pairs of orbits with nearby q values. Figure 
[16] depicts several pairs of orbits; each orbit in the center 
column may be approximated with the one in the left 
column. 



B. Comparisons Between RN and Schwarzschild 

The taxonomy also gives us a way to visually inspect 
orbits in different spacetimes to understand whether we 
can distinguish between them based solely on the dynam- 
ics of their periodic orbits. For example, an important 
question is whether it is possible to distinguish RN orbits 
from Schwarzschild orbits in this way. Figures |17| and [T8| 
depict periodic tables for the RN and Schwarzschild ge- 
omtries, respectively. By inspection, it is evident that 
high w orbits occur at higher energies in RN spacetime 
than their Schwarzschild counterparts. As a result, the 
apastra for RN orbits with w > are consistently larger 
than for Schwarzschild orbits. Comparing orbits in Fig- 
ure [18] to those in Figure [19| instead allows us to compare 
orbits with equal q for the same L; note however that 
here, the RN geometry is not extremal, as Q = 0.4. Also 
note that while the q = 1/2 orbit does not exist in the 
L = 0.4 Schwarzschild geometry, it can be found in the 
L = 3.8, Q = 0.4 spacetime. We still find that high w 
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FIG. 13: Periodic orbits with z = 1, 2, 3, and 4 for L = 3.8, Q = 0.5, 
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FIG. 14: The (w,V,z) = (w,0, 1) orbits for L = 3.5, Q = 0.5. From left to right, these are (0,0, 1), (1,0, 1), and 
(2,0, 1). Orbits with w > 3 are indistinguishable at this scale because the additional "whirls" are densely packed 
in r. 
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FIG. 15: Two z = 4, w = 1 orbits with different v values. 
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FIG. 16: A periodic table of orbits; in each column the upper orbit may be 
approximated by the corresponding orbit in the second row. 
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orbits occur at higher energies in RN spacetime and that, 
again, the apastra for these orbits are consistently larger 
than those of their Schwarzschild counterparts. 

We may use energy level diagrams, such as the one in 
Figure [20| to understand the relationship between q and 
E in each spacetime (for a general discussion, see [10j). 
These diagrams make it clear that RN orbits do appear at 
a higher set of energies. Next we want to look at how the 
spacing between these energies varies. Figure [21] shows 
the differences between values of E for successive orbits 
in the RN and Schwarzschild cases. 

We find that not only are the RN energies collectively 
higher, they spacing between each pair of energies is also 
consistently larger than in the Schwarzschild case. Fur- 
thermore, the magnitude of difference between spacings 
in each geometry increases. The ratio of the energy dif- 
ference in the RN and Schwarzschild q = 2 and q = 2g 
orbits is 3.52, but this rises to 23.41 when we compare 
the difference between the q = 6 and q — 6^ orbits in 
each geometry. 



VI. CONCLUSIONS 

We have applied the taxonomy of [9] to the Reissner- 
Nordstrom spacetime and demonstrated that within cer- 
tain bounds for Q and L, the orbital dynamics can be 
characterized using zoom-whirl behavior. These bounds 
were calculated for both the charged and uncharged par- 
ticle cases. Furthermore it was demonstrated that the 
regions in which we find bounded orbits in RN spacetime 
are analogous to those in Schwarzschild spacetime, and 
that we do not find multiple stable circular orbits for any 
choice of parameters within our predefined bounds. Ap- 
plying the taxonomy to RN spacetime for both charged 
and uncharged particles enables us to differentiate be- 
tween the set of periodic orbits in each geometry based 
on their orbital dynamics. We find not only that RN 
orbits occur at higher energies than their Schwarzschild 
counterparts, but that this is a behavior that persists for 
various Q. Furthermore, we find that RN orbits of a given 
q are more eccentric than their Schwarzschild equivalents. 



Appendix A: Hamiltonian formulation of RN 
geodesic motion 

For completeness we present a Hamiltonian formula- 
tion of RN motion. These were the equations integrated 
to generate the orbits in the paper. We begin with 
the Lagrangian density for a free particle in Reissner- 
Nordstrom spacetime: 



C 



a nP 



^gapq q 



(Al) 



The q a are dimensionless coordinates, as we have as- 
sumed that the orbiting particle is of unit mass. Fur- 
thermore, for timelike trajectories, C — —1/2. Using 



Pa 



dC 
dq~* 



we obtain the components of the momentum: 

Pt = lH = - Ai 
dC „ 

Pr = — = A 



at 

dC 
Or 
d£ 

dd 

oc 



-if 



■2a 



2 • 2 , 

r sin I 



The Hamiltonian is then defined 



n = P ^r - £ 



(A2) 

(A3) 
(A4) 
(A5) 
(A6) 

(A7) 



Like the q a , the particle's 4- momentum is also dimen- 
sionless and is here equivalent to the 4-velocity (since /z, 
the particle mass, is set equal to 1). So we may write 
Equation (A7) as 



n = -g ap Pa p P . 



If we compute % we find that 
1 



H 



A- 1 (-At) 2 + A- 1 (Ar) 2 



2/ 2a\2 i -2 • -2 at 2 • 2 al\2 

• " - n + r sin (r sin do) 



(A8) 



(A9) 
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Because each quantity in the Hamiltonian is just half the 
contraction of the 4-momentum, it is identical to the La- 
grangian, which is not surprising because our Lagrangian 
and Hamiltonian contain only kinetic terms [5]. 

Next we wish to plug our Hamiltonian into Hamilton's 
equations, 



% 
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dp % 



Pi 
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dqi' 



(A10) 
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FIG. 17: All extremal RN z — 1, 2, 3 orbits with w = for the first column, w — 1 for the middle 
column, and w = 2 for the last column. For all orbits, L = 3.3 and Q — 1.0. Orbits increase in energy 
from top to bottom and left to right. Note that the first and second entries in the first column are 
blank because the q — + j and q = + | orbits are inaccessible. 
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All Schwarzschild z — 1, 2, 3 orbits with w — for the first column, w — 1 for the middle 
and w = 2 for the last column. Orbits increase in energy from top to bottom and left to 
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FIG. 19: All RN z = 1, 2, 3 orbits with w = for the first column, w — 1 for the middle column, and 
w — 2 for the last column. Orbits increase in energy from top to bottom and left to right. 
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FIG. 20: Energy levels for Schwarzschild orbits (first row) with L = 3.8 and RN Orbits (second row) with L — 3.3, Q — f.0, 
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FIG. 21: Differences in energies at which we find orbits of 
various g in the RN (marked with "+") and Schwarzschild 
(marked with "•") spacetimes. 
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which requires that we rewrite H in terms of the pi. 
Hamilton's equations give us explicitly 
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The equations above were used in the numerical results 
presented in the body of the paper. 
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